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Abstract 

We present a detailed spectral and timing analysis of the High Mass X-ray Binary (HMXB) 4U 1909+07 
with INTEGRAL and RXTE . 4U 1909+07 is a persistent accreting X-ray pulsar with a period of approxi- 
mately 605 s. The period changes erratically consistent with a random walk expected for a wind accreting 
system. INTEGRAL detects the source with an average of 2.4 cps (corresponding to 15mCrab), but 
sometimes exhibits flaring activity up to 50 cps (i.e. 300mCrab). The strongly energy dependent pulse 
profile shows a double peaked structure at low energies and only a single narrow peak at energies above 
20 keV. The phase averaged spectrum is well described by a powerlaw modified at higher energies by 
an exponential cutoff and photoelectric absorption at low energies. In addition at 6.4 keV a strong iron 
fluorescence line and at lower energies a black body component are present. We performed phase re- 
solved spectroscopy to study the pulse phase dependence of the spectral parameters: while most spectral 
parameters are constant within uncertainties, the blackbody normalization and the cutoff folding energy 
vary strongly with phase. 
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1 Introduction 

4U 1909+07 (also known as X1908+07) was discovered with the Uhuru satellite as 3U 1912+07 (in the 
third Uhuru catalog, Giacconi et al. 1974). Since then, the source was detected with most X-ray instruments. 
The nature of the system was, however, not clear for 30 years. Only in 2000, Wen et al. (2000) found a 
stable 4.4 d period in RXTE ASM data, which was interpreted as the orbital period of a binary system. The 
photoelectric absorption was found to vary by a factor of at least 3 with orbital phase (Levine et al. 2004). 
Such a behavior can be well described by a spherical wind model and an inclination of 54° < i < 70°, 
depending on the parameters of the wind model. Furthermore, 4U 1909+07 shows irrgular fhring activity 
due to inhomogeneities in the wind of the donor star, which lead to different accretion rates. 

In 2005, Morel & Grosdidier (2005) detected an OB star in the near infrared at the location of the X-ray 
source, thus confi rming that the system is indeed a High Mass X-ray Binary (HMXB). They estimated the 
distance to the source to 7 kpc (Morel & Grosdidier 2005). 

Using a pointed RXTE - PCA observation. Levine et al. (2004) discovered a period of ~605 s, explained 
as the pulse period of a rotating neutron star with an offset magnetic axis. Levine et al. (2004) were thus 
able to refi ne the binary orbit parameters. They obtained A4 = 9-31 M e and R* < 22 R 0 for the mass and 
radius of the companion star assuming the canonical mass of 1 .4 M G for the neutron star. 

In this paper, we present the data and data reduction methods in Sect. 2. In Sect. 3 we analyze the pulse 
period evolution and the pulse profi les and in Sect. 4 we address the spectral properties of the source. We 
summarize and discuss our results in Sect. 5. 


2 Data and Observations 

We used all publicly available data of 4U 1909+07 from the International Gamma Ray Laboratory (INTE- 
GRAL , Winkler et al. 2003) and the Rossi X-ray Timing Explorer (RXTE, Jahoda et al. 1 996). 

INTEGRAL has three X-ray science instruments, which provide coverage from 3keV up to lOMeV: 
the imager IBIS/ISGRI (Ubertini et al. 2003, 20keV to 800 keV) with moderate energy resolution and a 
large effective area, the spectrometer SPI (20keV to lOMeV; Vedrenne et al. 2003) with excellent spectral 
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Figure 1: Closeup on an ISGRI lightcurve of 4U 1909+07 in the 20-40 keV with 20 s time resolution. The 
red line shows the folded pulse profi le. 

resolution suitable for the analysis of nuclear lines , and the X-ray monitor JEM-X (Lund et al . 2003 , 3 keV to 
35 keV) . Thanks to ISGRI’s large fi eld of view of almost 30“ x 30° more than 6 Msec of data on 4U 1 909+07 
exist, mostly from the INTEGRAL core programme (Winkler 2001) and the GRS 1915+105 monitoring 
(Rodriguez & Bodaghee 2008), however,4U 1909+07 was almost never in the fi eld of view of JEM-X, as 
no observations were pointed directly at 4U 1 909+07 . 

Among the detected sources in the fi eld of view are prominent sources like GRS 1915+105 or SS 433. 
In total, twelve sources with a significance > 6<r were detected in the field of view, with 4U 1909+07 
being the fourth brightest source. In addition to the standard pipelines of the Offline Scientific Analysis 
software (OSA) version 7.0 for images and spectra, the iiJLight distributed with OSA was used to obtain 
light curves with a higher temporal resolution. We used a time resolution of 20 s to maintain an acceptable 
signal-to-noise ratio and to obtain a good enough time resolution to measure the pulse period. 

In addition to the INTEGRAL observations, we also used the observations of 4U 1909+07 performed 
by RXTE. Both major instruments of RXTE, the Proportional Counter Array (PCA; Jahoda et al. 2006), 
sensitive in the 2-60keV energy band and the High Energy X-Ray Timing Experiment (HEXTE; Rothschild 
etal. 1998), sensitive between 1 5-250 keV, were fully operational at that time. Both instruments havealarge 
effective area and a high temporal resolution, but no imaging capabilities. In total, 196ksec of RXTE data 
on 4U 1909+07 are available. 


3 Timing Analysis 

3.1 Pulse period evolution 

Apart from the discovery of the pulsations based on RXTE data in 2000 and 2004 (Levine et al. 2004), no 
other measurements of the pulse period of 4U 1909+07 are available. The extensive archival INTEGRAL 
data are therefore optimal to track the evolution of the pulse period since 2003. We split the lightcurve 
into segments between 300ksec and 800ksec length to determine the period accurately. We then used 
epoch folding (Leahy et al. 1983) to determine the pulse period for each of these segments individually. 
4U 1909+07 shows very strong pulse-to-pulse variations such that the actual pulse can not be seen by the 
naked eye (see, e.g.. Fig. 1). Only after folding 500-1000 pulses, a stable pulse profi le emerges which 
allows for a reliable pulse period determination. The pulse profi le itself, however, is remarkably stable on 
longer time scales. A similar effect is seen in Vela X-l , although the overall luminosity is much higher in 
that source (Staubert et al. 1980). 

The evolution of the pulse period since 2001 is shown in Fig. 2. The overall behvior is best described 
by a random walk model , but also a spin-up trend between March 2003 and April 2006 and a spin-down 
trend from April 2006 to October 2007 could describe the data. To confi rm the random walk behavior, 
we implemented the algorithm proposed by de Kool & Anzer (1993): this algorithm evaluates the relative 
change in the pulse period P over different time intervals St between individual measurements. In the case 
of a perfect random walk, the result will be a straight line with a slope of 0.5 in the logd&i - log* space, 
where to is the angular velocity: to = 2n/P. Fig. 3 shows the corresponding diagram for 4U 1909+07. 
Superimposed is a line with slope of 0.5. shifted in y-direction to fi t the data. It is obvious that the line 
matches the data very well. The uncertainties in Fig. 3 include uncertainties of the determined pulsed 
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Figure 2: Evolution of the pulse period of 4U 1909+07. The historic RXTE data points of Levine et al. 
(2004) are shown as red crosses, our measurements obtained with INTEGRAL are shown as stars. 



Figure 3 : Pulse period evolution in log Sto - log 5t space , as proposed by de Kool & Anzer ( 1 993) . 


periods and also the uneven and coarse sampling. 

3.2 Pulsprofiles 

As shown by Levine et al. (2004), the 3.7 - 17 keV pulse profi le shows two distinct peaks, with the second 
peak being slightly broader than the fi rst peak and having a complex shape with two subpeaks at phase 
0 = 0.85 and <j> = 0. 1 . To study the energy dependence of the pulse profi le, we extracted pulse profi les from 
PCA data in different energy bands, using a period of P = 604.685 s and 32 phase bins. Thanks to the large 
effective area of the PCA instrument, high quality pulse profi les could be extracted in 30 narrow energy 
bands. Three energy bands are shown as an example in the upper three panels of Fig. 4. The pulse profi les 
show a smooth transition from a two peaked profi le at low energies to a single peak profi le at high energies. 
At energies below 5 keV the secondary peak is broader and stronger than the primary peak (Fig. 4a). With 
increasing energy both peaks become at fi rst equally strong and more clearly separated and then the relative 
strength of the secondary peak declines further (Fig. 4b ,c) . To obtain the pulse profile at high energies 
between 20keV and 40 keV, we used INTEGRAL data taken in 2004. We used the same epoch as for the 
PCA analysis, but a period of P = 604.747 s, as determined by our analysis of the INTEGRAL -data. In this 
energy band, the secondary peak does not exist anymore while the primary pulse is clearly seen and has a 
sharp peak (Fig. 4d). The deep minimum around phase 0.3, however, is not energy dependent. 


4 Spectral analysis 

4.1 Phase averaged spectrum 

The pulse phase averaged spectrum of 4U 1909+07 is best described by a powerlaw continuum attenuated 
by photoelectric absorption at low energies and an exponential turnover at high energies similar to most 




O *o 

c 

« 

2 - 


.a) 3.3-4. Ike V 

______ 

- 

■'-+ J 

Li l* 3 

_b) 9.1-10.3keV 


c) 1 4.9-1 6.2keV 


Id) 20~40keV * 


IS 

1 \ \ 


Figure 4: Energy resolved pulse profi les with RXTE PCA (a - c) and INTEGRAL ISGRI (d). The profi les 
are shown twice for clarity. Note that the RXTE and INTEGRAL profi les are not strictly phase aligned. 


accreting X-ray pulsars. Levine et al. (2004) modeled the spectrum using bremsstrahlung which describes 
the data equally well. The turnover at high energies is often modeled with the cutoffpl,highecut,or the 
fdcut (Tanaka 1986). The npex model (Mihara 1995) is more complex as it involves a second powerlaw. 
Furthermore at 6.4 keV an iron fhorescence line is present. 

We discarded all observations between orbital phase 0.88 < d > OI b < 0.12, as the Ah is dramatically 
increased during this part of the orbit (Levine et al. 2004). We applied the typical continum models (see 
above) to the RXTE and INTEGRAL data (see Table 1). All models can describe the data almost equally 
well. In the case of the fdcut model, however, the cutoff energy is set to <1 keV, thus effectively removing 
the cutoff. We therefore did not use the fdcut model any more for 4U 1 909+07 . Independent of the applied 
continuum models, however, a soft excess below 10 keV is evident, which can be very well modeled using a 
blackbody with a temperature of kT « 1 .4 keV. Fig. 5 shows the spectrum and the best fi t cutof fpl model. 

4.2 Phase resolved spectra 

Since the INTEGRAL ISGRI data did not provide high enough statistics for high resolution phase resolved 
spectroscopy, we only used the RXTE data. We divided the PCA and HEXTE data into 7 phasebins (see 
Fig. 6a). We applied a cutoffpl with an additional black body component to all phase bins. Neither the 
photoelectric absorption nor the power law index F varies with pulse phase. The cutoff folding energy, 
however, is highest in the primary peak with ~ 19keV, while it is only ~14keV during the secondary peak, 
explaining why the primary peak is much stronger at higher energies. 

The black body component is strongly variable: the normalization changes by a factor of ~3. The 
blackbody is strongest between the two peaks and weakest in the rise and maximum of the primary peak. 
The temperature kT of the black body on the other hand does not change with pulse phase. 


5 Discussion & Conclusion 

We have presented a detailed study with INTEGRAL and RXTE of 4U 1909+07. We have shown that 
the pulse period changes strongly with time. The evolution of the period is best described by a random 
walk. Such a behavior has been seen in many other HMXBs like Vela X-l and is a strong indicator for a 
direct accretion from the stellar wind of the optical companion without a persistent accretion disk (Ghosh 

6 Lamb 1979). An accretion disk would provide a more constant transfer of angular momentum and 
thus a long-term spin-up or spin-down trend, as seen in other sources like 4U 1907+09 (Fritz et al. 2006). 
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Table 1 : Fit parameters for various model of the phase averaged spectrum of 4U 1909+07 . 


Model 

parameter 

Cutoffpl 

Cutoffpl 

±bbody 

Highecut 

A red 

Ni{ [atoms cnr 2 l 

r 1 

Fe cr 

Fe Energy [keV] 
Fe Norm 

1.76 
15 3 +0 ' 6 
1 63 +o13 

fit? 2 

64^ 

U -7-0.2 

(1-0 x 10 

1.01 

4 7 + 1.6 
+• '- 1.9 

0'96-ao§ 

0-28!;® 

6 40 +d o4 
- 0.06 

‘ 3 (0.57^jj) x 10~ 3 

1.07 
4 8 +0 - 9 

1 ' 37 + 0-03 

-Q.Q8 
n 41+0.07 
U -^ 1 -Q ,06 
6 39+0.04 
- 0.03 

(0.78!«;» 7 ) x 10 -3 

Model 

parameter 

Highecut 

+bbody 

NPEX 

NPEX 

+bbody 

Nn [atoms cm 2 ] 

r 2 

Fe cr 

Fe Energy [keV) 
Fe Norm 

0.91 

6 -8!a8 
1.32 ±0.10 
0.2 ±0.2 

6 39+0 04 
- 0.05 

(51 3 ) x 10- 4 

1.44 

1 ,7 +0 ' 7 

0.36!«;?|/-3.08^ 

0 4J+0.07 

637^:85 

(0.80:°° 7 )xl0” 3 

0.93 

5+1 

0.80^“ /-2.0 3 
n 2Y+ 0 - 12 
6.40 ± 0.05 
(0.6!®; 3 )x 10' 3 



Figure 5: a Combined RXTE PCA (red) and HEXTE (blue) and INTEGRAL ISGRI (green) spectrum. The 
best fi t cutoffpl-model is also shown, b Residuals of the best fi t showing no evident deviation. 
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Phase 


Figure 6: Spectral parameters of phase-resolved spectra b Pulse profi le in the 1 3.2-14.5 keV energy range. 
The different shaded areas indicate the phasebins . b Photoelectric absorption column /v’h . c Power 1 aw index 
I', d Folding energy, e Black body normalization, f Black body temperature, and g width of the Iron line. 
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Since the accreting neutron star is strongly magnetized (as indicated by strong coherent pulsations), the 
accreted matter has to couple to the magnetic fi eld thus inhibiting the formation of an accretion disk close 
to the neutron star. The required B fi eld strength is of the order of ltf 2 G. Such magnetic fi eld strengths are 
commonly observed in neutron stars, especially pulsars. In such a strong magnetic fi eld, however, cyclotron 
resonant scattering features (CRSFs) are expected to be present in the X-ray spectrum, as in the spectrum 
of Vela X-l (Kreykenbohm et al. 2002) or MXB 0656-072 (McBride et al. 2006). Unfortunately, we found 
no evidence for a CRSF in 4U 1909+07 and are thus unable to determine the strength of the magnetic fi eld. 
However, the absence of CRSFs in the X-ray spectrum does not rule out the presence of a strong magnetic 
field in 4U 1909+07, as CRSFs can be filled up again by photon spawning and depend strongly on the 
geometry of the accretion column and on the viewing angle (Schonherr et al. 2007). 

The two distinct peaks in the pulse profi le at low energies suggest that the neutron star is accreting 
on both magnetic poles. However, the dramatic change of the pulse shape with energy suggests that the 
physical conditions in the accretion columns above the two poles might be different. The exact geometry 
of the accretion column, however, depends strongly on the way the matter couples to the magnetic fi eld and 
is highly uncertain (Basko & Sunyaev 1976; Meszaros 1984). Besides a simple fi lied funnel other possible 
confi gurations exist including a hollow or a partly hollow funnel and several small funnels. 

These models, however, have to be taken with a pinch of salt, as the behavior of relativisitc plasma in 
magnetic fi elds with B w 1(J 2 G is still only little understood. Additionally gravitational light bending must 
be taken into account when analyzing pulse profi les and their origin (Kraus et al. 2003). 
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